TITLE: 

TYPE: 

NUMBER OF WORDS: 

TEMPORARY STORAGE: 

DURATION: 

DESCRIPTION: 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



AUXILIARY 

LIBRARY ROUTINE E 8 - 307 



Integration by Weddle's Rule 

Closed 

30 

0, 1 

3.4 millisec. 

This routine uses Weddle's Rule to approximate an integral: 



bsf f(x) to ~ h [f + 5f l + f 2 + 6f 3 + h + f5 5 + f 6 ] 



where 



f f (a + nh) 
n x 



ENTRY: 



b-a 
is the value of f (x) at x = a + nh where h = -g- and is 

tabulated at memory location t + n. There are seven function 

values in all. 

This routine is entered by means of the instructions 
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26" ~ 
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where t is the location of the first tabulated function 
value f Q . 

When control is returned to the right side of p + 1, the 
approximation to the integral will be in both the 
accumulator and quotient registers. 



- 2 - 



ACCURACY: 



The main error terra in Weddle's rule is 



EXAMPLE: 



NOTE: 



. 1 h T f (6) 

150 



(J) ; 



so that results will be exact for a polynomial of fifth 
degree or less. 

Use of this routine to compute 
.1 







| sin jtx dx = - = .318309886 



gave the result 

.31830127. 

Using the same ordinates, integration by Simpson's Rule 

(code E - 2) gave a value of .31844727, with about 15 

times the error of this routine. 

1. Because the coefficients in Weddle's Rule are small 
integers, the intermediate terms in the calculation can 
be generated without using multiplication or division 

with the result that this routine is at least three and one- 
half times faster than the Simpson Method used in code E-2 for 
the same number of ordinates. The main advantage of this 
routine is thus its speed. 

2. The limits a and b are not needed, since the result 
is independent of a linear change of variable. 

3. Because of the method used, accuracy may be lost for 
very small integrands . 



DATE July 5, i960 



PROGRAMMED BY 



ohn. Ehrman 



APPROVED B Y , ^\^}h^A.ycU^V 



tr 



LOCATION 


ORDER 


NOTES 


PAGE 1 E 8 




OOK 









K5F 
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Plant link address 
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1»6 8L 








lA 22L 


Plant f Q address 
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46 17L 








IA 22L 


Plant f.. address 
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k6 12L 








10 20F 


Plant f 2 address 
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L4 8L 








42 23L 


Plant f- address 




5 


L4 8L 


«-) 






42 9L 


Plant f^ address 
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L4 8L 
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Plant fV address 




7 


L4 8L 
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Plant f> address 
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L5 F 








10 IF 
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40 F 
L5 F 
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10 IF 
L4 F 
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10 IF 


^^ 






40 IF 


iA jf + * J 
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L5 F 
10 IF 
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kO F 
L5 F 
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10 IF 
Lk F 
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10 IF 
L4 IF 
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1/16 Jf Q + f 2 + f^ + 
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L5 F 






10 IF 
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40 F 
L5 F 
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10 IF 
Lk F 
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10 IF 
40 F 
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Lk F 


5/16 & + fl 
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Lk IF 


L <J 






10 IF 


1/32 £f Q + 5^ + f 2 + + f k + 5f 5 + 
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40 IF 
L5 F 
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10 IF 

kO F 
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10 IF 
L4 F 


3A f 3 
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Lk IF 


6/32 f 3 




27 


JO 20L 
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clear Q (roughly), divide by 5/8 
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22 F 


exit via link 
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constant = 5/8 





